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Abstract. - We introduce a hard-core lattice-gas model on generalized Bethe lattices and 
investigate analytically and numerically its compaction behavior. If compactified slowly, the 
system undergoes a first-order crystallization transition. If compactified much faster, the system 
stays in a meta-stable liquid state and undergoes a glass transition under further compaction. 
We show that this behavior is induced by geometrical frustration which appears due to the 
existence of short loops in the generalized Bethe lattices. We also compare our results to 
numerical simulations of a three-dimensional analog of the model. 



The last decades have seen a great interest in understanding the structural glass transition 
^12]. The latter is given by a dramatic dynamical slowing down which prevents the systems 
from reaching the equilibrium state at experimental time scales. For glass- forming liquids, e.g., 
the glass transition is usually identified by the temperature where the viscosity exceeds 10 13 
poise. A central role in theoretical approaches is played by fully connected spin-glass models 
with multi-spin interactions 012] ■ These models in fact show a purely dynamical transition 
towards a glassy state followed by a thermodynamic transition at lower temperature. The 
dynamical transition, even if not strictly present in real materials due to activated processes, 
is believed to present some similarities to the structural glass transition. 

The analogies between fully connected spin glasses and structural glasses are, however, of 
purely phenomcnological nature, missing a microscopic verification. The main differences for 
spin-glass models are, that they have multi-spin interaction, which in addition are disordered 
(quenched disorder) and long-ranged (fully connected models). So, in particular, these models 
do not have any crystalline state. 

Recently, some of these points could be cured by considering models based on Bethe lattices 
or Husimi trees, which are the analog of Bethe lattices for multi-spin interactions: On finite 
length scales, these lattices have no inhomogeneities. The corresponding models undergo a 
crystallization transition which was found to be of first order in the most interesting cases. By 
cooling or compactifying these models relatively fast, the system stays in the meta-stable liquid 
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(or paramagnetic) phase. In analogy to their fully-connected counterparts, they undergo a 
dynamical glass transition at even lower temperature, which is followed by a thermodynamic 
transition 0J|S]. The latter cannot be observed in simulations (or experiments) since the 
system does not equilibrate below the dynamical transition. These models still have, however, 
the problem of having multi-spin interactions (which, in some cases, are introduced by local 
constraints concerning more than two particles, as seen in 0). Similar models were also 
discussed in the context of compactificd granular matter [5]. 

Here we present a lattice-gas model with pure two-particle interactions. Defined on a usual, 
i.e., locally tree-like Bethe lattice, the model does neither show a first-order crystallization 
nor a glass transition. If, however, the Bethe lattice is generalized to include short loops, a 
dynamical behavior reminiscent of the glass transition in liquids of hard particles or Lennard 
Jones liquids [7] is found. This glass transition is induced by geometrical frustration 8 which 
prohibits constructing a globally dense configuration from purely local rules. 

As the original Bethe lattice, also the generalized Bethe lattices studied here are of mean- 
field nature. In order to compare our results to more realistic models, we introduce a three- 
dimensional lattice gas showing the same kind of geometrical frustration. We study this model 
be means of Monte-Carlo (MC) simulations: The behavior of the model shares many features 
with the mean-field case, even if they are in general much less pronounced. 




Fig. 1 Fig. 2 



Fig. 1 - Part of the generalized Bethe lattice. The unit cell is equivalent to a clique of p vertices 
(here p = 4). Each site is contained in K + 1 unit cells (here K = 3). 

Fig. 2 - Density p as a function of the chemical potential p, for p = 4, K = 3. The full lines give 
analytical results for the liquid and the crystalline phases. The spinodal and the crystallization points, 
as well as the dynamic and static glass transition are marked by vertical lines (from left to right). 
Results are compared with numerical compaction curves for random generalized Bethe lattices of size 
N — 999, averaged over 100 graphs. Here Sp = 0.2 was used. Inset: Compaction rate dependence for 
high values of p for 71mc ranging from 5000 (top) to 10 (bottom). 

The hard-sphere lattice-gas model can be defined on any graph G = (V,E). Each site 
i e y can be occupied by at most one particle, xi £ {0, 1}. Particles are assumed to have a 
hard core of radius one, i.e. neighboring sites can never be occupied simultaneously: XiXj = 
if {i,j} £ E. Note that this defines a two-particle interaction. If the graph is a tree, i.e. if 
it does not contain any loops, dense particle packings can be generated from a simple local 



Martin Weigt and Alexander K. Hartmann : Glassy lattice GAS 



3 



rule: Starting with one particle, all neighboring sites have to be empty due to the hard- 
core constraint. Particles can again be put to all second neighbors and so on. Geometrical 
frustration enters via loops. If, e.g., some of the second neighbors of a site are directly 
connected by an edge, they cannot be occupied simultaneously. The resulting arbitrarity 
in selecting occupied second neighbors finally leads to globally disordered packings of sub- 
maximal density. 

To investigate the influence of geometrical frustration, let us consider the following gener- 
alized Bethe lattices: They are constructed from basic units consisting of p-cliques, i.e. fully 
connected subgraphs of p vertices. From these basic units, a graph is constructed by joining 
K + 1 cliques in every vertex, see Fig^ A usual Bethe lattice is obtained for p = 2 where 
the basic cliques are just simple edges, and the graph is locally tree-like. For p > 2, the graph 
contains triangles, i.e. local loops, giving rise to geometrical frustration. 

The thermodynamic properties of the model are given by the grand-canonical partition 
function 

5= E e ^ iXi II (!) 

{xi}G{0,l} N {i,j}£E 

The last product includes the hard-core constraint, it equals one for allowed packings, and 
zero if any two neighboring sites are occupied simultaneously. The chemical potential fi is 
coupled to the total particle number and can be used to tune the particle density. In our 
system, we study the compaction of the model which can be achieved by increasing /i. In 
particular, densest packings are reached in the limit fi — ► oo. 

The model can be solved most easily using the cavity method which, in the liquid as 
well as in the crystalline state, reduces to the Bethe-Peierls iterative method. Note, however, 
that the latter is modified with respect to the presence of short loops. On a coarse-grained 
level, where every local clique is considered as one object, the (hyper-)graph becomes locally 
tree-like and therefore solvable by an iterative procedure. In the following, we use this coarse- 
grained tree-like structure and speak about trees, branches etc. only in this sense. The cliques 
are numbered by v = 1, (K+l)N/p. We closely follow the presentation and notation of 
cf. also [S]. 

Let us consider a branch of the graph which is rooted in vertex i. One of the cliques 
containing vertex i is not part of the branch, we denote it by v. We first calculate the 
partition function S /i(i, v) of this branch where i is fixed to be empty (xi = 0) or occupied 
[xi = 1). Let us denote the other K cliques containing i by v\, vk-, and the further vertices 
contained in clique Vk by ji,fc, jp-i,k, k = 1, K. Due to the hard-sphere constraint, every 
clique is allowed to carry at most one particle, we thus find 



A' 



fc=i 



K 



p— 1 p— 1 

J\ 3o0'o,fcj v k) + E S i(^,fc> v k) Y[ E o(jb,k, Vk) 
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Introducing cavity fields [ih^ — ln(Si(i, i>)/So(i, v)), these iterate according to 



K 



*:=i 




(3) 
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i.e. they result from a linear superposition of the local chemical potential p and fields prop- 
agated along the cliques; and they take values h\ G [1 — K — ln(p — 1 + e~ M )/^, 1]. The 
true effective field acting on a vertex is given by the contributions propagated from all neigh- 
bors, i.e. from all K + 1 cliques containing the vertex, H j = 1 + u (^ji k > k ) e 
[-K — ln(p — f + e~^)/p, 1]. Note that the r.h.s. of this equation still contains the fields hj 
coming from single branches only. The distribution of these fields serves as an order param- 
eter and allows to identify liquid, crystalline and glassy states. Furthermore, its knowledge 
allows to calculate the grand canonical potential with £ denoting the set of all p-cliques 
in the graph, — In 5 = — K J2i In (f + e^ ffi ) +XLe£ ^ n + Siei> e ^ hi ) > as weu as * ne avera g e 
particle density p per site, Np = Xi — ^ e ^ H< /(I + e ^ Hi ) ■ 

The liquid phase, which is expected to be globally stable for small chemical potential p, is 
characterized by its translational invariance. The effective fields h\ and consequently Hi do 
not depend on the site index i and are thus given by the homogeneous fixed point of Eq. 

For K > 2, this translational invariance breaks down in the crystalline phase: On each 
clique, there is exactly one site of high average density, all the other (p — f ) vertices have 
homogeneous lower density. There are thus two fields M 1 ' < h^ 2 \ which are a solution of 

= l + Ku(hW,hW,...,hW) 

= I + Ku(hS l \...M 1) ) ■ (4) 

For p = 2, i.e. usual Bethc lattices, the crystalline solution appears continuously due to a 
local instability of the liquid solution, the solidification is thus of second order, and the model 
easily crystallizes even if compactified relatively fast 

In the remaining part of the paper we therefore concentrate to the more interesting case 
p > 3, i.e. to graphs showing geometrical frustration due to short loops. There the crystalline 
solution appears discontinuously at a spinodal point p s . Still, its grand canonical potential 
is larger than the one of the liquid phase. The two potentials cross at the crystallization 
transition p c > p s , which thus describes a first-order transition being accompanied by a 
discontinuous density jump. Sending p — » oo, the crystalline system quickly approaches the 
maximal density 1/p. For the example p = 4, K = 3, which is presented in Fig. |2 the 
transition points are given by p s ~ 1.73 and p c ~ 2.18. 

The liquid phase stays, however, meta-stable even above p c describing thus a super-cooled 
liquid. If we compactify the system rather quickly, it initially stays liquid. However, at higher 
density, the system even falls out of this local equilibrium: An exponential number of glassy 
states of smaller density (and higher grand canonical potential) appears discontinuously at a 
dynamical transition p& and traps the system. At p ls \, a sub-exponential number of glassy 
states finally reaches lower grand canonical potential than the liquid one. In the equivalent 
replica approach, this corresponds to one-step replica symmetry breaking. For p = 4 and 
K = 3 we find pd ~ 5.7 and p rs b — 6.1, i.e. they are located far beyond the crystalliza- 
tion transition; at pa there are about exp(0. 00887V) metastable states, but the clustering 
phenomenon becomes more pronounced for larger p and K . 

Technically, this can be established using the cavity approach which generalizes the Bethe- 
Peierls solution to more than one state The main difference to the liquid solution is that, 
in different pure states, there are different cavity fields h\. They can be characterized by their 
histograms P"(K) which have to be determined self-consistently using the cavity equations, 
for technical details see the analogous case discussed in ^J. With Pi{h) being Dirac peaks 
in the liquid phase, the dynamical transition point p^ can be identified by the discontinuous 
appearance of a non-trivial solution for P"(h) which imidiately covers the full intervall of 
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allowed fields. The static transition /i rs b follows from the point where the replica-symmetry 
broken grand canonical potential becomes dominant with respect to the liquid one |12| . 

For comparison, we have performed numerical simulations obtained on fixed-connectivity 
random graphs built from 4-cliques. From an analytical point of view they share the liquid 
and the glassy state with the generalized Bethe lattices discussed above. Only the crystalline 
state is inconsistent with the random character of the graph entering via disordered loops of 
length 0(N), and it cannot be observed in the simulations. The simulation works as follows. 
For each MC step, a site i is selected randomly. With probability p — 0.5 a MOVE (M) step 
is performed, and with probability 1 — p an EXCHANGE (EX) step: 

M If site i is empty and has exactly one occupied neighbor, the particle is moved to i. In 
all other cases, the configuration remains unchanged. 

EX If the site is occupied, the particle is removed with probability exp(— fi). If the site is 
empty, and all neighboring sites are empty, a particle is placed on i. 

Note that in this way detailed balance is fulfilled, and the equilibrium distribution is the one 
corresponding to Eq. (Q. All simulations start with fi = 0. Then /i is increased in steps of 
Sp = 0.1. For every value of the chemical potential, we perform umcN Monte-Carlo steps. 

In Fig. [21 the analytic findings are compared with numerical results. Note that in principle 
the value of /i<j can be inferred from the point where the numerical curve deviates from the 
liquid one in the limit of infinitesimal slow compaction (%c — > °°)- Even if our data turned 
out to be incompatible with every tested form of fits (as a function of tt-mc); this point appears 
to be consistent with the analytical prediction which is marked by a vertical line in the figure. 

As already mentioned, also the generalized Bethe lattices are or mean-field nature. Whereas 
a similar glassy behavior was found also numerically in a two-dimensional model of cross- 
shaped particles ^2], it is useful to test our theoretical predictions at least qualitatively in 
a finite-dimensional model which is closer to our mean-field model. In particular it should 
be characterized by the same kind of geometrical frustration. We thus consider a three- 
dimensional lattice with N = L 3 sites. It is generated from a tetrahedron, i.e. a 4-clique, 
which is translated along its edges. Again each site is contained in K + 1 = 4 tctrahedra. 
We impose periodic boundary conditions. On top of this lattice, we again consider our lattice 
gas of hard spheres of radius one, and we numerically study its behavior by means of MC 
simulations up to size N = 100 3 . 

As in its mean-field counter-part, the system is found to be in a liquid phase for small /i: 
All particles are free to move around, they are not localized, and the time-averaged single- 
site density is homogeneous. For slow compaction, we observe a discontinuous crystallization 
transition at // ~ 2, cf. Fig. At this point, the symmetry between sites breaks down, 
every tetrahedron has one site of higher than average density, and three of lower density. The 
average density exhibits a discontinuity and quickly approaches 0.25, i.e. its theoretically 
maximal value, if we further increase /i. The first order character of this transition can be 
numerically tested by decreasing \i again. Doing this we observe a hysteresis, cf. Fig. 01 the 
system stays in the crystalline state even for fj, < 2. 

If, starting from the liquid state at low fj,, we increase /i quite fast, crystallization is avoided. 
The system stays at lower than the crystalline density. This phenomenon appears in analogy 
to the super-cooled, meta-stable liquid observed in the mean-field model. At higher fi, the 
system starts to freeze into a glassy state, being characterized by a density which is much 
smaller than 0.25. This asymptotic density depends heavily on the compaction rate. 

There are, however, some crucial differences between Figs. and |31 If we consider the 
compaction curves in the mean-field model, cf. Fig. they overlap even beyond the crystal- 
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Fig. 3 Fig. 4 

Fig. 3 - Density p of the lattice model as a function of the chemical potential p for timc = 
1,2,5,25, 1000 Monte-Carlo sweeps per value of the chemical potential (from bottom to top). All 
umc sweeps, the chemical potential is increased by 0.1. The data are for one run of a lattice with 
N — 50 3 sites. The solid line represents the analytical result for the crystalline phase on the gener- 
alized Bethe lattice which can be considered as a simple mean-field approximation. Inset: Density 
p as a function of p for slowly increasing (timc = 500) and subsequently fast decreasing (timc = 5) 
chemical potential. 

Fig. 4 - Density p of the lattice model as a function of the chemical potential p for umc = 2, when 
the chemical potential is first decreased every timc sweeps by 0.1 until p a t op is reached. Then the 
system continues to evolve at fixed value of p. The inset shows the density as a function of time for 
smaller p s t op - The arrows indicate the times from where the chemical potential is kept fixed. The 
horizontal lines denote the equilibrium densities of the crystalline state at ^i s to P . 

lization transition, until the (local) equilibration time start to be larger than the characteristic 
compaction time, and the system fall out of the meta-stable equilibrium. If we stop to increase 
p at some point below the glass transition, the density evolves towards the one of the liquid 
phase, i.e. the latter is a true metastable phase. This seems not to be true for the finite- 
dimensional lattice. The compaction rate dependence sets in directly at the crystallization 
transition, and there is no pronounced gap between a meta-stable liquid state and the globally 
stable solid one. In Fig. 0]we report the following experiment: The system is compactified 
quite fast by increasing the chemical potential up to some value p s top which is chosen beyond 
the crystallization point. Then the system is allowed to evolve at fixed p. As can be seen 
in the inset for a system of size N = 50 3 , the model still equilibrates towards the crystalline 
state (given by the horizontal lines), even if this equilibration time is much larger than the 
equilibration time in the liquid phase. In the main part of Fig. ^ we show the system size 
dependence of this effect: The equilibration time grows with N, and probably diverges in the 
thermodynamic limit. This growth is not exponential in N, i.e. there are no extensive energy 
barriers to be overcome as in the generalized Bethe-lattice model. 

To conclude, we have introduced a hard-core lattice-gas model with only two-particle 
interactions. On a generalized Bethe lattice, it exhibits a crystalline phase as well as glassy 
behavior if super-cooled sufficiently fast. In contrast to usual Bethe lattices, the generalized 
ones have short loops of odd length. These loops, and the resulting geometrical frustration, are 
responsible for the glassy behavior of the model which is similar to the behavior known from 
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models with multi-spin interactions: If super-cooled, the system falls out of equilibrium at a 
dynamical transition which is accompanied by the appearance of an exponential number of 
states. Also in our model, this transition is not of thermodynamic nature, a replica-symmetry 
breaking equilibrium transition occurs only at higher chemical potential. It is quite interesting 
that similar phenomena can also be found in a completely different model considered in |14j : 
There, geometrical frustration is introduced into a field theoretical setting in form of a small 
non-Abelian background gauge field. In a sense, these different approaches complement and 
re-enforce each other. 

Furthermore, we have studied numerically a three-dimensional variant of the model ex- 
hibiting the same local parameters, like sizes of the short loops and number of neighbors. The 
behavior of the lattice model is close to the mean-field type system except that the lattice 
variant does not have a true metastable phase which remains stable at the timescale of the 
simulation. 

Another conclusion concerns the typical numerical hardness of the vertex-cover problem, 
which is equivalent to finding dense hard-sphere packings 9 : Graphs containing small loops, 
or geometrical frustration, are expected to be much harder to cover with a small vertex cover 
size than locally tree-like graphs like random graphs. This opens the interesting question in 
how far local structures may influence the global hardness of a combinatorial problem. 
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